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SUMMARY 

The  continuum  low  density  hypersonic  flow  of  a  chemically 
react  mg  and  radiating  gas  over  a  blunt  body  is  analyzed  by  an 
approximate  technique  based  on  a  recent  alternative  formulation  of 
the  well-known  method  oi  integral  relations.  Specifically, 
approximations  are  introduced  for  the  streamwise  variations  of 
several  flow  properties,  while  detailed  radial  profiles,  (consistent 
with  the  stipulated  procedure)  are  sought  by  numerical  analysis 
along  normals  to  the  bod'y  (essentially  normals  to  streamlines); 
the  particular  formulation  allows  careful  analysis  of  the  large 
changes  in  physico-chemical  properties  ot  the  gas  between  body 
and  she  *  as  are  to  be  encountered  at  the  low  Reynolds  numbers 
where  viscous  effects  are  important  over  the  entire  shock  layer. 
Detailed  derivations  or  the  simplified  system  cc  ordinary  differ¬ 
entia!  equations  and  oi  the  oour.dary  conditions  governing  the 
problem  is  presented.  An  outline  of  numerical  procedures  for 
applications  13  also  given. 
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I.  INTRODUCTION 

The  problem  of  observables  associated  with  reentry  of  hyper¬ 
sonic  vehicles  has  received  considerable  attention  in  recent  years . 
This  problem  nas  strong  physico-chemical  as  well  as  fluid  mechan¬ 
ical  overtones.  Applications  require  detailed  and  accurate 
quantitative  predictions  of  flow  properties;  such  level  of  precision 
can  only  be  obtained  by  numerical  methods.  A  large  body  of  related 
literature  has  recently  been  set  forth;  however,  most  analyses  are 
pertinent  to  high  Reynolds  number  flows  where  distinct  shock, 
inviscid  and  boundary  layer  regions  can  be  recognized.  The  range 
of  applicability  of  those  analyses  is  limited  to  altitudes  below 
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approximately  220  kft.  for  typical  body  dimensions  (R=l  ft),  and 
hypersonic  reentry  conditions.  At  the  higher  altitudes  different 
flow  regimes  are  encountered,  which  involve  increasingly  complex 
analyses.  In  an  upward  sweep  of  the  trajectory,  a  continuum 
description  of  the  flow  can  be  applied  to  study  the  initial 
vorticity  interaction,  viscous  layer,  and  incipient  merged  layer 
regimes  while  the  range  of  lower  Reynolds  numbers  must  be  analyzed 
on  the  basis  of  either  transitional  or  free  molecule  flow  theory. 

The  present  study  is  concerned  with  the  first  three  (con¬ 
tinuum)  regimes  mentioned  above.  Even  within  these  limits  analysis 
becomes  very  cumbersome,  particularly  when  effects  of  chemical 
reactions  and  radiation  are  included.  Current  state  of  the  art 
permits  detailed  numerical  investigation  of  either  strictly 
inviscid  flows  (Ref.  2)  or  strictly  boundary  layer  flows  (Ref.  3) 
including  coupled  chemical  reactions  and,  in  preliminary  fashion, 
radiation*  (Ref.  4).  Low  density  situations  have  been  studied 
to  limited  extent  and  only  in  connection  with  flows  of  ideal 
gases;  although  related  analyses  rely  heavily  on  extensive  machine 

For  missile  and  lunar  vehicle  reentry,  the  effect  of  self 
absorption  can  become  significant  in  some  part  of  the  shock 
layer  (Ref.  4);  however,  in  this  preliminary  study  we  assume  the 
thin  gas  approximation  to  be  valid  over  the  entire  flow  field. 
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Calculations  results  are  only  available  for  the  stagnation  point 
on  a  sphere  (Ref.  5)  and  for  flows  about  blunt  bodies  within  the 
Newtonian  approximation  (Ref.  6).  The  limitations  of  those 
theories  are  well  known. 

The  present  analysis  was  undertaken  having  in  mind  the 
strict  requirements  for  prediction  of  observables,  definition  of 
the  range  of  possible  observation,  and  interpretation  thereof. 
However,  in  view  of  the  exploratory  nature  of  the  task  an  approx¬ 
imate  method  of  analysis  has  been  sought,  presumably  retaining 
the  essential  features  of  the  physical  situation. 

In  the  oresence  of  viscous  effects  and  of  chemical  reaction 
extending  to  the  entire  shock  '  dyer  careful  consideration  must 
be  given  to  distributions  of  state  properties  and  composition  in 
the  direction  normal  to  the  body  surface;  on  the  basis  of  this 
consideration,  we  have  adopted  an  integral  method  formulation 
based  on  approximate  description  of  the  distribution  of  flow 
properties  in  the  direction  parallel  to  the  body  surface 
(Refs .  7  and  8) . 

Details  of  the  investigation  are  presented  in  the  following 
sequence:  Section  II)  The  Governing  Partial  Differential 

Equations,  Boundary  Conditions  and  the  Approximations  Based  on 
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Orders  of  Magnitude  Analysis;  Section  III)  The  Integral 
Approximations  and  the  Resulting  Ordinary  Differential  Equations 
and  Boundary  Conditions;  Section  IV)  Outline  of  Numerical 
Procedure  Including  Iterative  Techniques  Successfully  Proven  in 
Connection  with  other  Problems. 

The  present  analysis  represents  the  first  phase  of  a 
research  program  on  reentry  observables  in  low  density  continuum 
flows  currently  under  way  at  GASL. 
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1 1  *  GOVERNING  EQUATIONS  AND  BOUNDARY  CONDITIONS 

The  equations  governing  the  two-dimensional  or  axisymmetric 
flow  of  a  viscous,  chemically  reacting  and  radiating  gas,  in  a 
body  oriented  coordinate  system  (x,yt$)  are: 


Continuity : 
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i  =  1,2.. . N-l 
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Where: 
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body  curvature 

0  two-dimensional  flow 
l  1  axisymmetric  flow 
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Fuithermore  for  Le  =  Pr  =  1  and  optically  thin  gas 
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The  system  of  4+N-l  equations  is  complemented  by  an  equation  of 
state  in  the  form: 


h  =  Mp.p.OL) 


(8) 


by  a  viscosity  law; 


M  = 


(9) 


and  by  the  condition: 


i=l 


(10) 
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The  corresponding  unknowns  in  the  original  system  are  the 
two  components  of  velocity  u,  v,  the  pressure  p,  the  density  p  and 
(N-l)  mass  fractions. 

Two  sets  of  boundary  conditions  are  imposed:  one  set  at  the 
body  and  one  at  the  "outer  edge"  of  the  shock  layer;  we  specify 
all  the  boundary  conditions  in  a  form  appropriate  for  the  "viscous 
layer  regime,"  namely:  1)  no  slip  and  temperature  jump  phenomena 
are  present  at  the  surface  of  the  body,  and  2)  the  shock  can  be 
treated  as  a  surface  of  discontinuity. 

The  boundary  conditions  at  the  body  are: 


%  *  vb  =  0 

h,  ‘  V*’ 


Surface 

catalyticity :  1)  if  catalytic 


in  equilibrium 


>(11) 


2)  if  not  catalytic  -♦impermeability  con¬ 
ditions  for  all 
species  except 
surface  material 

/ 


The  boundary  conditions  at  the  shock  are  obtained  from  the 
standard  Rankine-Hugoniot  relations  in  terms  of  free  stream  prop¬ 
erties, of  local  shock  incl inat ion,  and  of  a  prescribed  composition 
of  the  mixture  on  the  downstream  side;  formally 
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u  =  u 
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The  total  orde  *  of  the  original  system  is  [7+2(N-l)]?  the 
[ 7+2 (N-l ) ]  boundary  conditions  (11)  and  (12)  are  however  insuf¬ 
ficient  since  the  shock  geometry  is  also  unknown.  The  diffi¬ 
culty  is  removed  upon  examination  of  the  continuity  equation  at 
the  body?  this  provides  the  additional  boundary  conditions, 
namely 


(v  h  -  0  (11a) 

V  b 

The  original  system  of  Eqs.  ( 1 )  —  ( 5 )  may  be  simplified  con¬ 
siderably  for  specific  application  to  high  Mach  number  flows. 

Under  those  conditions  the  orders  of  magnitude  of  the  basic 


properties  in  the  shock  layer  are: 
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Accordingly  the  following  simplified  form  of  the  equations  is 
obtained  when  terms  of  order  higher  than  ( <R^  )  are  neglected 
(Ref.  5): 
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Momentum  y : 
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The  system  (14)  represents  the  basis  of  the  present  analysis. 
All  subsequent  developments  are  aimed  at  developing  a  r.ethod  of 
solution.  In  this  connection  we  begin  by  introducing  new 
independent  variables 

s  =  x 


n  =  y/fl 


(15) 
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with 


6  =  6 (s) 


(16) 


the  shock  detachment  distance;  the  domain  of  integration  is  thus 
transformed  into  a  rectangle. 

Performing  the  transformation 
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and  recasting  the  equations  in  divergence  form  (desirable  for 
implementation  of  the  integral  approach)  we  obtain: 
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Momentum  y: 


The  boundary  conditions  are  stiil  expressed  by  (11)  and  (12) 
intended  to  apply  at  n=0  and  at  n=l  respectively. 
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HI.  THE  INTEGRAL  APPROACH 

Metuodology  for  solution  of  the  simplified  system  of  partial 
differential  Eqs.  (18)  is  established  only  for  a  particular  formula 
ti°n  of  the  problem  (the  inverse  problem  where  the  shock  shape  is 
specified  and  the  associated  body  is  determined)  and  involves  con- 
si  lerable  labor.  The  practical  problem  of  analyzing  the  flow 
about  a  given  configuration  can  be  solved  within  practical  limits 
of  numerical  effort  only  in  approximate  fashion.  A  well  known 
procedure  to  this  effect  is  represented  by  the  integral  method  in 
its  various  formulations  (Refs.  7  and  8).  in  essence  this  method 
reduces  the  system  of  governing  partial  differential  equations  to 
an  approximating  system  of  ordinary  differential  equations  by 
introducing  assumed  distributions  of  flow  properties  along  one 
c  ordinafe  direction.  in  the  present  report  we  choose  to  make 
assumptions  about  the  variation  of  flow  properties  in  the  stream- 
wise  direction  since,  in  the  presence  of  viscous  effects  and  of 
chemical  reactions  extending  over  the  entire  shock  layer,  we  are 
particularly  interested  in  assessing  the  effect  of  transport 
phenomena  on  cross-stream  distributions  of  physico-chemical  prop¬ 
erties  of  the  gas. 

We  begin  implementation  of  the  aforenoted  integral  approach 
by  integrating  Eqs.  (18)  between  two  general  stations  sj  and  ss 
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to  obtain  the  integro-dif ferential  equations 
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Reduction  of  the  (3+N)  Eqs.  (19)  to  ordinary  differential  form 
follows  immediately  upon  stipulation  of  the  s-distribution  of 
certain  flow  quantities.  We  point  out  that  only  (3+N)  quantities 
can  be  subject  of  independent  description;  of  these,  two  must  be 
velocities,  two  state  properties,  and  (N-l)  species.  We  assume 
polynomial  laws  for  the  following: 
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pu  =  Ai  s  +  A3  s3  +  •  •  • 
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where  the  coefficients  A,  B,  C,  D,  E  are  functions  of  n.  The 
selected  quantities  are  all  equal  to  zero  at  the  axis  and  the 
assumed  s  dependencies  satisfy  the  symmetry  conditions;  there¬ 
fore,  the  number  of  terms  in  each  polynomial  is  equal  to  the  number 
of  control  stations  away  from  the  axis.  Consistently  with  this 
procedure  the  properties  at  the  axis  are  found  from  ratios  of  the 
laws  (20)  and  not  from  the  independent  solution  of  the  stagnation 
point  partial  differential  equations  (Ref.  5);  for  example: 
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and  so  on  for  all  other  properties.  We  also  notice  that  in  the 
case  of  one  control  station  and  one  term  in  the  laws  (20),  the 
proposed  procedure  is  equivalent  to  assuming  flow  similarity. 
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In  addition  to  the  four  independent  quantities  (20)  other 
accessory  properties  are  cast  in  a  polynomial  form  to  permit  a  more 
convenient  numerical  integration  of  the  final  system;  naturally  the 
degree  of  these  polynomials  is  determined  by  the  constitutive  rela¬ 
tions  to  be  satisfied  by  the  quantities  on  hand  [e.g.  Eqs.  (6)] 
consistent  with  the  postulated  distributions  (20).  Thus,  we  set 
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consistent  with  symmetry  requirements  and  specific  definitions. 
Upon  substitution  of  (20)  into  (19)  and  (22)  one  obtains  t 
simultaneous  nonlinear  ordinary  differential  equations  of  the 
first  order  which  can  be  cast  in  the  form 
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where 


CT  X  =  1,2,3. ..t 


t  is  defined  in  terms  of  the  number  z  of  control  stations 

(s  =  const.  =  si  ,  sa  .  .  .  )  ,  the  number  of  chemical  species  N  and 

the  magnitude  of  the  highest  powers  nu  appearing  in  (22) 

t  =  z ( 3+N)  +  m!  +  ma  +  m3  +  (N-l)m*  (23a) 


P  (n)  denotes  the  general  coefficient  (A.,B  ,  etc.)  in  the 
A  11 

stipulated  distributions  (20)  and  (22),  a^  and  are  known 

functions  of  s  ,  n  and  P.  to  be  determined  case  by  case  in  accord 

A 

with  the  number  z  of  stations  considered.  Of  the  t  Eqs.  (23), 
z(3+N)  follow  from  the  conservation  laws  (19),  while  the  remaining 
[mj  +ma +m3  +  (N-IJra* 1  are  obtained  by  setting  to  zero  the  coefficients 
of  all  powers  of  s  in  (22).  The  t  boundary  conditions  for  the 
system  (23)  may  readily  be  obtained  in  the  following  way:  z(3+n) 
conditions  are  obtained  by  simultaneous  consideration  of  (11), 

(12)  and  (20);  [mi +mg +m3 + (N-l ) m* ]  additional  conditions  are 
derived  therefrcm  by  substitution  into  (22)  and  subsequent 
setting  to  zero  of  the  coefficients  of  all  powers  of  s.  Recog¬ 
nition  of  additional  z  unknowns,  namely  the  z  coefficients  in  the 
law  describing  the  shock  shape,  does  not  alter  the  order  of  the 
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system  and  merely  requires  implementation  of  the  conditions  (11a) 
at  the  z  control  stations. 

In  summary,  an  integral  method  analysis  of  the  viscous  flow 
of  a  mixture  of  N  reactant  gases  over  either  two-dimensional  or 
axisymmetric  bodies,  by  recognizing  z  regions  (strips)  within  the 
domain  of  interest,  involves  the  solution  of  a  system  of  t 
ordinary  differential  equations  of  first  order  (23)  plus  z  alge¬ 
braic  equations,  subject  to  t  boundary  conditions  [t  defined  by 
(23a)].  The  detailed  derivation  of  expressions  for  coefficients 
and  forcing  functions  in  the  system  must  be  carried  out  case  by 
case  (depending  upon  the  selected  number  z  and  the  thermodynamic 
behavior  of  the  mixture);  however,  some  rules  of  internal  con¬ 
sistency  must  be  respected  in  that  process  as  described  below. 

To  be  consistent  with  the  stipulation  (20)  various  flow  prop¬ 
erties  of  interest  must  be  described  by  the  relations: 


pv 


pv 


a 


(puv)  (pu) 
(PU*  ) 


l£uvJl 

(Pu5) 


(pua. ) 

(pu) 


(24a) 

(24b) 


a 

i 


(24c) 
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h 


=  Ifiuhi 

(Pu) 


(24d) 


H 


=  h  +  f  (us+va)  =  h  +  \  {^yi-  Mpus)  +  (pv3)] 


(24e) 


P 


phv [s.n) 


(puh) (pu) 
(pu) 


ff(s,n) 


( 24  f ) 


The  formal  expression  for  the  pressure  p  is  suggested  by  consid¬ 
eration  of  ideal  gas  flows  wherein  if  =  const;  in  the  more  general 
case  one  has 

tr  (s,  n)  =  ir(p, h.fltj  <25 


As  an  example  of  consistent  derivation  let  us  examine  a  compli 
cated  term  such  as  that  involving  the  pressure  in  the  x-momentum 
Eq.  (19b); 


ds 


(26) 


The  consistent  law  for  the  shock  is 


6  60  +  6i  sa  +  6a  s*  + 


Thus 


rj 

6 


<36  „  ,  v 

"  R°(s’ 


+  jn  Ri  (s) 


(27) 


(28) 
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with 


Ro  (s) 


1 

6 


d6 

ds 


/cos  0„ 


Ri  (s)  = 


dfi 

ds 


cos  0, 


The  appropriate  sequence  of  transformations  of  (26)  then  proceeds 
along  the  following  lines: 


i_  rSs  r£  d6  3  rSar 

Sn  J  6  ds  P  ds  “  3n  Ro(s)+jn  Rj  (s)  x 

si  s*  J 


ff(s,n)ds 


s*  *r  (s ,  n)  ds  | 


3n 


Ro (s)+jn  Ri (s) 

Di  +B4S® 


T?  (s,  n)ds 
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S(DaAi  fAflD,  ) 
dn 


R0  (s) +jn  Ri  (s) 
B3+B4sa 


s3  it  (s,n)dp 


+  txl  rSs 

dn  J 

si 


Rq (s)+jn  R^  (s) 
Bg  +B4  S*5 


s4  rr  (s ,  n)  ds 
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(Ba+B*s^r 


ff(s,n) ds 
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s4  tr(s,n)ds 
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All.  the  integrals  in  (30)  except  the  last  one  may  be  evaluated  in 
straightforward  fashion  to  obtain  associated  contributions  to  the 
coefficients  a  .  ;  the  last  integral  is  reduced  to  the  desired 

aX 

form  by  means  of  the  additional  relationship: 


p  =  R0 (s)+jn  Ri (s) 


3a  Aj  +Di  )  S3  +Da  An  S* 
Bj  +B4  s3 


(32) 


All  coefficients  and  forcing  functions  in  (23)  may  con¬ 
sistently  be  determined  by  algebraic  transformations  as  detailed 
at  (26)  through  (32)  above. 


26 


Solution  of  the  system  (23)  may  only  be  sought  by  numerical 
means;  proper  formulation  thereof  requires  a  preliminary  assess¬ 
ment  of  the  general  features  of  the  integral  curves  to  be  obtained. 
In  this  connection  we  observe  that  the  integral  formulation  adopted 
here  does  not  exhibit  critical  points  of  the  type  present  in  the 
standard  integral  method;  hence  the  integration  is  straightforward. 
The  absence  of  singularities  follows  from  the  a  priori  assumed 
regular  behavior  at  the  sonic  velocity  as  manifested  either  by 
the  prescribed  pressure  distribution  on  the  body  or  by  the  pre¬ 
scribed  shock  shape  (see  the  alternative  procedures  in  Section  IV) . 
The  singularity  in  the  standard  formulation  represents  a  mathe¬ 
matical  manifestation  of  the  continuity  requirement;  in  that  case 
we  start  by  assuming  a  shock  detachment  distance  60  at  the  axis 
and  compute  increment  thereof.  Here  we  do  not  say  anything  about 
6;  it  is  completely  free  to  satisfy  the  continuity  requirement 
consistent  with  the  assumed  regular  behavior. 

Consideration  of  viscous  effects,  specifically  the  no-slip 
condition  at  the  wall,  in  the  present  formulation  renders  the 
equations  singular  at  n=0;  however,  this  singularity  is  regular 
and  integration  in  its  neighborhood  may  readily  be  performed 
along  the  lines  detailed  in  Ref.  5. 
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outline  for  numerical  procedure 

In  analogy  with  the  inviscid  blunt  body  problem  two  possi¬ 
bilities  are  open  for  the  problem  on  hand:  the  direct  and  the 
inverse  method. 

Analysis  by  the  direct  method  proceeds  along  the  following 
lines:  a)  the  body  geometry  and  the  body  temperature  are  given? 
b)  (N+3)  distributions  of  flow  properties  (polynomials  in  s)  arc 
estimated  at  n=0;  these  are,  for  example,  the  shock  geometry  6(s), 
the  pressure  p,  the  shear  stress  [whereby  a  relation  between 

!  ^  jn  obtained  ,  and  either  the  composition  if  the 

wall  is  not  catalytic,  or  the  normal  derivative  of  the  N  species 

°  r  1  l  pu  In  ^  t^lC  wail  is  catalytic;  c)  the  remaining 
inputs  for  the  integration  are  consistently  obtained  from  the 
study  of  the  equations  in  the  neighborhood  of  the  body?  specifi¬ 
cally,  one  requires  v^O  to  satisfy  continuity  and  one  determines 
pn  to  satisfy  y  momentum  equation  consistent  with  the  pre¬ 
scribed  p(s)  ?  d)  numerical  integration  of  (23)  is  performed  from 
n=0  to  n=l?  o)  compliance  of  the  results  of  integration  at  n=l 

w.’ th  the  boundary  (jump)  conditions  for  p,  p,  u,  v  and  a  is 

i 


•k 

The  aforenoted  set  of  inputs  is  equivalent  to  prescribing  the 
heat  transfer  since  h  =h  (p  ,n  a  ) 

n  n'*n  yn‘  in'  * 
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inspected?  and  f)  if  the  latter  conditions  are  not  satisfied 
initial  estimated  values  are  modified  and  an  iterative  solution 
performed  following  criteria  such  as  are  discussed  later  in  this 
section . 

Analysis  by  the  inverse  method  proceeds  along  the  following 

lines:  a)  the  shock  geometry  is  given;  b)  (4+N)  distributions 

of  flow  properties  (e.g.  6,u  ,p  ,h  ,0t  )  are  estimated  and  pre- 

^  n  n  n  in 

scribed  at  the  outer  edge  of  the  shock  layer  (n=l)  and  initial 
values  for  integration  established  accordingly;  c)  integration 
is  carried  out  from  n=l  to  n= 0;  d)  compliance  with  (4+N)  condi¬ 
tions  on  the  body  (namely  u~v=v^-0,  h=h^  and  either  a  or 
consistent  witn  surface  cataly ticity)  is  inspected;  and  e)  itera¬ 
tion  is  pursued  if  not  all  conditions  are  satisfied  at  n=0. 

It  must  be  concluded  that  either  method  exhibits  comparable 
numerical  difficulties  for  the  problem  on  hand;  indeed  a  large 
number  of  quantities  [either  (3+N)  or  (4+N)  must  be  estimated 
and  iterated  upon  for  either  direct  or  inverse  approach. 

The  required  iterations  may  be  carried  out  -.n  organized 
fashion  by  assuming  linearization  of  differences  between  solutions 
characterized  by  change  in  one  input  parameter  only.  To  be 
specific,  a  first  run  is  performed  correspondxng  to  an  initial  set 
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of  (3+N)  estimated  quantities;  subsequently  (3+N)  additional  runs 
are  performed  each  characterized  by  a  set  of  initial  conditions 
differing  from  the  first  one  for  a  single  quantity.  In  the 
hypothesis  that  the  error  can  be  linearized  a  system  of  linear 
algebraic  equations  in  the  unknown  A  [ j=l, 2, • • • (3+N) ]  correc¬ 
tions  to  be  imposed  on  the  estimated  initial  conditions  may  be 
constructed;  specifically  one  equates  the  errors  observed  in  the 
first  run  to  a  weighted  linear  combination  of  the  errors  intro¬ 
duced  by  each  perturbation  to  obtain 


P. 

1 


3±N  dP . 


(33) 


th 

where  P  indicates  the  known  initial  error  property  for  the  i  ‘property 
and  the  partial  derivatives  are  obtained  by  differences  between 
the  results  of  each  "perturbation"  run  and  the  initial  one. 

Obviously  the  procedure  can  be  repeated  until  satisfactory  con¬ 
vergence  is  achieved. 

A  second  alternative  to  the  numerical  iteration  may  be 
sought  by  extension  to  the  problem  on  hand  of  Weil's  method  for 
solving  a  class  of  ordinary  differential  equations  usually 
encountered  in  boundary  layer  problems.  The  essence  of  that 


30 


method  is  the  approximate  solution  of  the  equations  by  iteration, 
satisfying  exactly  at  each  iteration  the  known  boundary  condi¬ 
tions.  Although  the  characteristics  of  convergence  can  not  be 
predicted,  an  approach  along  this  line  is  here  suggested  in  view 
of  its  simplicity,  its  adaptability  to  automatic  operation  and  its 
capability  to  provide  results  to  any  degree  of  accuracy.  The 
salient  procedural  points  are  clearly  illustrated  by  inspection 
of  the  following  model  solution  of  a  system  of  two  non-linear 
ordinary  differential  equations  of  the  first  order;  extension  to 
the  higher  order  system  (25)  is  immediate.  Consider  the  equations 


an 


dx 


+  a13 


dx 


Fi 


aai 


dll 

dx 


+  a33 


df? 

dx 


F3 


(34) 


with  allf  a13  ,  a3i,  a33  ,  Fi  ,  Fa  =  known  functions  of  fi  ,  f*  ,  x, 
subject  to  the  boundary  conditions 


at 


x  =  0  fi  =  fio 


(35) 


at 


x  =  1 


fa  —  fgl 
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Solve  the  system  (34)  for  the  highest  order  derivatives  (in  this 


case  first  order) : 


P.  -  f“-  F, 

_ aa  a 

ai,+a»,  f*- 

aaa 


(36) 


At  the  general  ntn  iteration  the  properties  (fx  )^  ^  and  (^a)n  ^ 
are  known  (in  the  initial  stage  estimated  profiles  satisfying  the 
known  boundary  conditions  are  used  as  inputs.  Write  Eq.  (36)  in 
the  form: 
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Upon  integration 


Cfi  (x)  Jn  =  fio  + 


Pi 
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The  boundary  conditions  at  x=l  define 


(  f  2  °  )  f  2  1  ~ 
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and,  therefore,  Eqs. 

(38) 

can  b 

files.  The  procedure 
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[(f„)n-(f3.>n-i] 
(  f  a  o  )  n 


.  The 


procedure  is  amenable  to  further  refinements:  in  particular 
when  the  solution  oscillates  an  over-relaxed  iteration  or  extrap¬ 
olation  can  be  used  to  expedite  convergence. 
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V.  CONCLUSIONS 

A  method  has  been  presented  for  the  study  of  the  continuum 
viscous  shock  layer  about  re-entry  bodies  as  it  is  encountered 
at  altitudes  approximately  between  300  and  200  kft  for  vehicles 
having  typical  dimensions  of  the  order  of  feet.  An  integral 
approach  has  been  used  which  permits  detailed  description  of 
distributions  of  flow  properties  along  normals  to  the  body 
surface  and,  thereby,  should  lead  to  reasonable  predictions  of 
observables.  The  present  report  has  been  concerned  with  the 
initial  phase  of  the  investigation,  namely  the  development  of  the 
method  of  analysis.  Applications  will  require  extensive  numerical 
work;  guide  lines  therefor  have  been  discussed. 
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Attn:  Mr.  E.J.  Little/AMSMI-RNR  (34-8010-42) 

USAMICOM 

Redstone  Arsenal 

Huntsville,  Alabama 


Institute  of  Science  &  Technology  (34-8010-44) 

The  University  of  Michigan 
Box  618 

Ann  Arbor,  Michigan 
Attn:  BAMIRAC  Library 

Raytheon  Company  (104-8010-45) 

Attn:  Mr.  A.  Truszanowski 

Sudbury,  Massachusetts 

R.C.A.  Laboratories  (34-8010-47) 

Moorestown,  New  Jersey 
Attn:  Mr.  Archie  Gold 


Advanced  Research  Projects  Agency 
Washington,  D.C.  20301 
Attn;  Dr.  B.  Fisher 
Dr.  C.D.  Fitz 
Mr.  R.  Zirkind 

Technical  Information  Office 
Mr,  C.  McLain,  Rm.  2B-263 


(34-8010-48) 

(34-8010-49) 

(34-8010-50) 

(34-8010-51) 

(34-8010-52) 


Foreign  Technology  Division  (64-8010-54) 

Wr ight-Patter son  Air  Force  Base,  Ohio 
Attn:  TD-FCS  (Schmidt) 


Conductron  Corporation  (64-8010-55) 

393  S.  Main  Street 
Ann  Arbor,  Michigan 
Attn:  Dr.  K.  Siegel 


1 

1 


1 


1 


1 


1 


1 

1 

1 

3 

1 

1 


1 
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Syracuse  University  Research  Corp.  (64-8010-56) 

Syracuse,  New  York 
Attn:  Paul  Howells 

Picatinny  Arsenal  (64-8010-57) 

U.  S.  Army 

Dover,  New  Jersey  07801 
Attn:  SMUPA-VL 

General  Motors  Defense  Research  Labs.  (64-8010-58) 
Aerospace  Operations  Department 
Santa  Barbara,  California 
Attn:  Dr.  Colin  Maiden 

Aeroelastic  and  Structures  Laboratory  (64-8010-59) 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 
Attn:  Dr.  Emmitt  A.  Witmer 

National  Engineering  Science  Company  (64-8010-60) 
711  South  Fair  Oaks  Avenue 
Pasadena,  California 
Attn:  Dr.  Richard  Tracey 

Douglas  Aircraft  Company  (84-8010-61) 

Missile  and  Space  Systems  Division 

Reentry  Systems  Advanced  Programs 

Santa  Monica,  California 

Attn:  V.  E.  Crosley 

TRW/Space  Technology  Laboratories  (64-8010-62) 

One  Space  Park 

Redondo  Beach,  California 

Attn:  Dr.  William  S.  Carlson 

Sperry  Gyroscope  Co.  (84-8010-63) 

Great  Neck,  New  Yolk 

Attn:  Mr.  K.  B.  Erkenbrack 


of  Cys. 

1 


1 


1 


1 


1 


1 


1 


1 
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Flui  Dyne  Engineering  Corp. 

5900  Olson  Memorial  Highway 
Minneapolis,  Minnesota  55422 
Attn:  Dr.  James  S.  Holdhusen 

Aerospace  Corporation 
P.0.  Box  1308 
San  Bernardino,  California 
Attn:  Mr.  K.  Pierce 

Mr.  G.  Milburn 

Deputy  for  Studies  and  Analyses 
Research  and  Technology  Division 
Attn:  SESOM  (Lt.  Col.  Greer) 

Wright-Patterson  AFB ,  Ohio 


No.  of  Cys. 

( 84-8010-64)  1 


(84-8010-65)  2 

(84-8010-66)  2 

(84-8010-67)  1 


UNCLASSIFIED 


UNCLASSIFIED 


